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We investigate the linear mechanical response of an aqueous foam, and its relation to the micro- 
scopic rearrangement dynamics of the bubble-packing structure. At rest, even though the foam is 
coarsening, the rheology is demonstrated to be linear. Under flow, shear-induced rearrangements 
compete with coarsening-induced rearrangements. The macroscopic consequences are captured by 
a novel rheological method in which a step-strain is superposed on an otherwise steady flow. 

PACS numbers: 82.70.Rr, 83.60.Rs, 83.80.Iz, 83.85.Cg 



Aqueous foams are tightly packed collections of gas 
bubbles separated by a continuous liquid phase P, |3| ■ 
Like elastic solids, bulk foams resist shear, completely 
unlike the gases and liquids from which the they are com- 
prised. The origin of this striking behavior is that the 
bubbles are jammed, unable to flow around one another 
and explore configuration space under thermal energy. 
Thus bubbles distort, rather than rearrange, when sub- 
jected to small shear deformations. The resulting extra 
internal gas-liquid surface area costs energy in proportion 
to the surface tension, and this provides a restoring force. 
The shear modulus is roughly surface tension divided by 
bubble radius, depending precisely on the volume frac- 
tion, e, of the continuous liquid phase ^ As a 
foam is made wetter, the bubbles become progressively 
rounder, and the shear modulus decreases. The elastic- 
ity completely vanishes at the point where the bubbles 
are close-packed spheres. This is one example of unjam- 
ming @. 

In this paper we explore other ways to unjam the bub- 
bles in a foam. Each could correspond to a different tra- 
jectory in a global jamming phase diagram 0. For ex- 
ample, one could imagine raising the temperature so that 
ksT is greater than cri?^, which would allow the bubbles 
to rearrange thermally like Brownian particles. However, 
a typical value is aB? /ks ~ 10^^ K; therefore, foams are 
athermal far-from-equilibrium systems and raising the 
temperature is not feasible. But there are at least two 
other ways to drive the system so that bubble rearrange- 
ments occur. One is through coarsening, the diffusion 
of gas from smaller to larger bubbles. As this proceeds, 
local stress inhomogeneities repeatedly build up to some 
theshold and relax by sudden avalanche- like local rear- 
rangements. Such microscopic dynamics have been ob- 
servedjDreviously via diffusing-wave spectroscopy (DWS) 
0, 1^ ||l£|. This raises a string of interesting questions. 
How do localized coarsening-induced rearrangements af- 
fect the linear shear rheology of bulk foams? Is the effect 
similar to thermally excited dynamic heterogeneities in 
glassy systems? An entirely different way to induce re- 
arrangments is by shear. Here local stresses also accrue 
and relax, but in a more correlated manner. These rear- 
rangment dynamics have similarly been studied by DWS 



|llL Il2l Il3| . How does the elastic character of foam van- 
ish as the shear rate is increased? Can we think of the 
two driving mechanisms as "heating" and ultimately un- 
jamming the sample? To make progress, we study the 
linear mechanical response of an aqueous foam at very 
long times and low frequencies. In a novel twist, we also 
study linear response during uniform shear flow. 

Our samples are a commercial aqueous foam, consist- 
ing of nearly-spherical polydisperse gas bubbles, 92% by 
volume, that are tightly packed in an aqueous solution 
of stearic acid and triethanolaminc (Gillette Company, 
Boston MA). Samples are surrounded by a water bath 
held at 25.0°C and are measured after the foam has aged 
for approximately 100 minutes. By this time, the aver- 
age bubble size is approximately 60 microns and is grow- 
ing reproducibly via coarsening 0|; stresses due to the 
loading process have also relaxed. Test durations are 
sufficiently short that gravitational liquid drainage and 
bubble coalescence are negligible. 

Our measurements are performed with a Paar Phys- 
ica UDS 200 rheometer, controlling the rotational speed 
and angular displacement of a solid cylinder whose axis 
is vertical and concentric with a fixed surrounding cup. 
The sample cell has an inner radius of 20.0 mm, and a 
4.1 mm gap; these dimensions ensure that the foam can 
be treated as a bulk material with uniform stress. To 
minimize end-flow effects, the 98 mm long inner cylinder 
is much shorter than the depth of the surrounding cup 
yet much longer than the gap width. Wall slip is pre- 
cluded by coating both the cylinder and cup with a fine- 
grade sandpaper. Samples are loaded through a 5 mm 
hole at the bottom of the cup, after lowering the cylin- 
der to the test position. DWS measurements indicate 
the absence of shear-banding and other secondary flows 
|l3| . To re-confirm, we compared with a similarly coated 
cone-and-plate cell, with a 10 cm radius and 10° cone 
angle. Identical rheology results were obtained for both 
cell configurations, implying that the imposed shear de- 
formation is uniform. Diffuse light transmission indicates 
that bubbles do not burst, and voids do not form, even 
under high shear 0, . 

Effect of coarsening on mechanical response. To quan- 
tify linear response, we measure both the complex shear 
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FIG. 1: Dynamic shear moduli of a coarsening foam, in both 
the (a) frequency and (b) time domains; symbols in (a) de- 
note the imposed strain amplitude. The dashed curves in 
(a-b) are a fit to G*(aj) — Go(l + y^iuj/ujn) and G{t) — 
Go(l + l/^/nuJni), predicted for nonaffine bubble motion in 
the absence of any time evolution. The solid curve in (b) is 
an empirical fit to G{t); appropriately transformed, it gives 
the solid curves mathing the storage and loss moduli in (a). 



modulus, G*(w) = G"(w) + iCJuAM, [T^, and the 
stress relaxation modulus, G{t) [l^llal- Our data are 
displayed in Fig. ^ Note that G'{uj) > G"{uj), mean- 
ing that the response is primarily elastic rather than 
dissipative. The frequency range for our G*{uj) data 
spans almost six decades, from 2t: divided by sample 
age (an absolute minimum below which measurement is 
not possible) up to a maximum set by limitations of the 
rheometer. Note that different amplitude strains give the 
same result, demonstrating absence of wall-slip and other 
geometry-dependent artifacts and hinting at linearity of 
response. A more stringent test of linearity is comparison 
with G{t), which should be related to G*{llj) by Fourier 
transform 0, 0|. The time range for our G{t) data 
spans over five decades, from the time needed to achieve 
the step strain up to the time beyond which stress is zero 
to within instrumental limits. An empirical fit to the 
G{t) data is shown in Fig. ^b) by a solid curve; this fit 
is transformed and plotted in Fig. ^a) over a frequency 
range corresponding to the time range of the fit. The 
agreement is very good, demonstrating conclusively that 
the sample is linear. This is further supported by an em- 
pirical fit to G*(w) data at high-w and the comparison of 
its transform with G{t) data at short-<. 

Let us now consider the frequency and time depen- 



dence of the moduli in Fig. ^ The fit for w > 5 rad/s 
is to the form G*(w) — Go(l + \/ iLu/ujn) with free pa- 
rameters Go = 2300 dyne/cm^ and = 156 rad/s. The 
former represents the static shear modulus, roughly sur- 
face tension divided by bubble size jj, U, ||5| . The lat- 
ter represents the effect of nonafRne deformation of the 
bubbles under shear due to local packing configurations 
that are strong or weak with respect to the shear direc- 
tion This fit, including the parameters, is consistent 
with the G*{uj) data in Ref. [ijj, where the frequency 
range (0.3 — 20 rad/s) was too small to fully demonstrate 
the functional form. According to the theory of 0| , the 
characteristic frequency is LUn oc Go/??oo, where the very- 
high frequency response is G*{uj) = irjaoUj. The numer- 
ical prefactor was not predicted; experimentally, it was 
found to depend on s and was not of order 1. To com- 
pare, the value is w„ « 600 rad/s for a e = 0.38 emulsion 
of 0.5 /im oil droplets in water [tH.!!^ (NB: by our defi- 
nition, ujn may be easily read off a plot by locating where 
G'(2a;„) = 2Go and G "(2u;„) = Go). 

The fit to Go(l + ^iujjun) fails for uj <h rad/s; the 
corresponding transform fails for i > 20 s. At longer 
times, the G(t) data decay slowly below Go, almost log- 
arithmically, over a few decades before relaxing more 
rapidly at around 1000 s. The transform of this final 
decay corresponds to the peak in G"{ijj) at 10"'^ rad/s. 
At lower frequencies, G*{lj) is unmeasurable; but since 
the integral of G{t) over all time is finite, the very-low 
frequency behavior is formally G"{uj) oc w as required 
by causality Thus the full frequency-dependence 

of G*{uj) for our foam is truly known and well-behaved. 
This resolves a long-standing puzzle raised by ear- 
lier measurements Il8l [Hiil where G*{uj) was roughly 
constant down to the lowest measured frequencies. 

All that remains is to understand the origin of the low- 
ijj j long-f behavior. We contend that evolution of the 
foam structure by coarsening is responsible. One clue 
is that the onset of deviation from the high-w fits cor- 
responds to the time Toq = 20 s given by DWS for the 
time between coarsening-induced rearrangements at each 
site. Another clue is that the final decay of G(t) and, 
equivalently, the peak in G"{lS) correspond to the sam- 
ple age. Since coarsening gives power-law growth, it takes 
of order the sample age for the structure to completely 
change. It is interesting that coarsening-induced rear- 
rangments relax microscopic coarsening-induced stress 
inhomogeneities far more quickly than the relaxation of 
macroscopically-imposed stress. Rather, the cumulative 
effect of many rearrangements and a change in bubble 
size is needed to relax global stress. The net result is a 
rheology that obeys linear response. This is remarkable 
given that the microscopic relaxation mechanism isn't 
thermal motion, but rather evolution. In effect, coars- 
ening unjams the foam, so that at low frequencies the 
rheology is G'[ui) <S G" (ui) = rjcu like an ordinary equi- 
librium liquid. 
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FIG. 2: (a) Stress relaxation moduli for foams sheared at 
various rates, as listed. This is given from the stress a{t) fol- 
lowing the superposition of a step-strain A7 on steady shear: 
G(t, 7) — [cr(t) — ct(0)]/A7. (b) The corresponding relaxation 
spectra, H{t,^) ~ dG{t,^)/d\nt. 



Effect of shear on mechanical response. Next we in- 
vestigate how the apphcation of shear causes a simi- 
lar loss of elasticity. For this we superimpose a small- 
amplitude step-strain, A7, on top of otherwise steady 
shear at rate 7. The resulting temporary increase in 
stress defines the transient stress relaxation modulus, 
G{t,j) = [a{t) — a{Q)]/Aj. A frequency-domain ver- 
sion of this technique was developed for polymers [2^. 
Example data are shown in Fig. |2Ia) . Note that as the 
strain rate increases, the elastic character of the foam 
melts away. There are two signs of this. First, the pro- 
gressively smaller intercept, G(0,7), means less elastic 
energy storage. Second, the progressively shorter decay 
time means more liquid-like dissipation. At high enough 
strain rates, where there is no transient storage and only 
dissipation, the foam behaves like an equilibrium liquid. 
Similar behavior has now been observed in simulations 
of a model foam . 

To further quantify the unjamming behavior appar- 
ent in Fig. [2Ja) , we first deduce the relaxation spectrum, 
H{t,j) « dG{t,j)/dlnt Results are shown in 

Fig. |2Jb). For low and zero strain rates, there are two 
competing relaxation processes, reflected by a broad peak 
at very late times and long tail of short-time modes of 
nearly equal weight. The former reflects the coarsening 
process, and the latter the nonafhne bubble motion. As 
the shear rate increases, the coarsening peak gradually 
falls and another peak gradually rises over the nonafhne 
tail at short times. Presumably this is due to shear- 



induced rearrangements. 

The salient features of the stress relaxation are shown 
in Fig. 121 vs strain rate. The first plot is of elastic storage, 
G'(0,7), the value when the superimposed step-strain is 
achieved (below about 0.1s). This decreases very slowly, 
and is nearly constant, for strain rates less than about 
0.05/s; for higher strain rates it abruptly vanishes. The 
second plot is of stress relaxation times. One such mea- 
sure is te, when stress falls to 1/e of the initial value. An- 
other measure is tp, where 11(1,"/) reaches a global max- 
imum. At zero and very low strain rates, these times are 
different since there are two competing relaxation mech- 
anisms (evolution and shear). For strain rates higher 
than about 3 x 10~'''/s, the stress relaxation is essentially 
exponential and the two relaxation times are hence indis- 
tinguishable. In this regime, shear completely dominates 
the relaxation. It is puzzling that the relaxation time 
decreases with increasing strain rate as 7"^/^, since on 
dimensional grounds one would have expected 

Now we may compare the macroscopic rheology with 
the nature of the microscopic bubble dynamics. Pre- 
viously we used DWS to measure the strain rate de- 
pendence of two microscopic time scales: Tq, the time 
between localized discrete rearrangements, and Ts, the 
time for adjacent scattering sites to convect apart by one 
wavelength of light The observed DWS data are 

reproduced by the dashed curves in Fig. Ol^b). For very 
low strain rates, below about 3 x 10~^/s, the rearrange- 
ments are discrete and the time between events, Tq, is 
not noticeably different from the quiescent value. The 
stress relaxation time is much much longer than Tq - in- 
dicating that the cooperative effect of many coarsening- 
induced rearrangements is required. For slightly higher 
rates, the time between events decreases, roughly as 
7""'^/^, just like the stress relaxation time. Since the 
stress relaxation time is now shorter than Tq, not ev- 
ery site in the foam must rearrange in order to relax 
the overall stress. Shear-induced events begin to domi- 
nate over coarsening-induced events at strain rates near 
7V = "ty/Toq = 0.025/s, where 7y = 0.05 is the yield strain 
[13. At still higher rates, events merge together and 
the flow becomes progressively more homogeneous and 
smooth. The crossover strain rate is 7^ = ly/rd = 0.5 /s, 
where = 0.1 s is the duration of rearrangements |l3| . 
Many physical quantities display a change in character 
above and below 7^ [l^. Indeed, just below 7^, the 
new rheological measures of elasticity vanish at precisely 
where we no longer can detect dicrete rearrangments (i.e. 
at the endpoint of the dashed curve for Tq). At higher 
strain rates, the bubble motion is dominated entirely by 
uniform shear, as characterized by the DWS time scale 
Ts ■ The onset of this correlated shearing motion of bub- 
bles (i.e. the starting point of the Ts curve) corresponds 
quite well to where G{0, 7) begins to drop. 

Finally, we note that the dramatic changes in the elas- 
tic character of the foam with strain rate are not reflected 
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FIG. 3; (a) Transient shear modulus, (b) transient shear re- 
laxation times, and (c) viscosity, all as a function of shear 
rate. The vertical dotted lines denote the characteristic 
shear rates set by the yield strain divided respectively by 
the time between rearrangements and the duration of rear- 
rangements in a quiescent foam. In (a), the the solid curve is 
Go exp[— 7/(0.1s~^)] where Go ~ 2300 dyne/cm^ is the shear 
modulus. In (b), the open diamonds indicate where H(t,^) is 
a global maximum and plusses indicate where G{t,'y) falls to 
1/e of its initial value; the solid line is a power law with ex- 
ponent of —1/2; the dashed curves indicate DWS time scales 
for rearrangements (tq) and shear (ts). 



very strongly in the viscosity of the foam. As seen in 
Fig. EJc), the viscosity decreases across the whole strain 
rate range, though not quite as fast as (which would 
have indicated constant stress) . There is at most a slight 
change in exponent at jm- Similar behavior is found in 
simulations [25j. Altogether this emphasizes the impor- 
tance of the superposition method to measure G(t, 7), 
since it provides a clear dramatic signature of the un- 
jamming transition. With this tool, and by comparison 
with DWS data, we have succeeded in connecting the 
nature of microscopic bubble dynamics with the result- 
ing macroscopic rheological behavior. We have thereby 
shown that the unjamming of foam can be accomplished 
both by time and by application of shear. The unjammed 
liquid-like state is very similar to what would be achieved 
by raising the temperature for a thermal system. An im- 
portant next step would be to deduce an effective temper- 



ature, recently shown in simulation of a sheared model 
foam to have many of the attributes of a true statistical 
mechanical temperature |26j . 
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